Symmetric multistep methods for constrained Hamiltonian systems
نویسندگان
چکیده
A method of choice for the long-time integration of constrained Hamiltonians systems is the Rattle algorithm. It is symmetric, symplectic, and nearly preserves the Hamiltonian, but it is only of order two and thus not efficient for high accuracy requirements. In this article we prove that certain symmetric linear multistep methods have the same qualitative behavior and can achieve an arbitrarily high order with a computational cost comparable to that of the Rattle algorithm. Mathematics Subject Classification (2010): 65L06, 65L80, 65P10
منابع مشابه
Symplectic and symmetric methods for the numerical solution of some mathematical models of celestial objects
In the last years, the theory of numerical methods for system of non-stiff and stiff ordinary differential equations has reached a certain maturity. So, there are many excellent codes which are based on Runge–Kutta methods, linear multistep methods, Obreshkov methods, hybrid methods or general linear methods. Although these methods have good accuracy and desirable stability properties such as A...
متن کاملReducing round-off errors in symmetric multistep methods
Certain symmetric linear multistep methods have an excellent long-time behavior when applied to second order Hamiltonian systems with or without constraints. For high accuracy computations round-off can be the dominating source of errors. This article shows how symmetric multistep methods should be implemented, so that round-off errors are minimized and propagate like a random walk.
متن کاملSymmetric schemes, time reversal symmetry and conservative methods for Hamiltonian systems(
It is important, when integrating numerically Hamiltonian problems, that the numerical methods retain some properties of the continuous problem such as the constants of motion and the time reversal symmetry. This may be a di cult task for multistep numerical methods. In the present paper we discuss the problem in the case of linear autonomous Hamiltonian systems and we show the equivalence amon...
متن کاملLong-Term Stability of Symmetric Partitioned Linear Multistep Methods
Long-time integration of Hamiltonian systems is an important issue in many applications – for example the planetary motion in astronomy or simulations in molecular dynamics. Symplectic and symmetric one-step methods are known to have favorable numerical features like near energy preservation over long times and at most linear error growth for nearly integrable systems. This work studies the sui...
متن کاملDilations, models, scattering and spectral problems of 1D discrete Hamiltonian systems
In this paper, the maximal dissipative extensions of a symmetric singular 1D discrete Hamiltonian operator with maximal deficiency indices (2,2) (in limit-circle cases at ±∞) and acting in the Hilbert space ℓ_{Ω}²(Z;C²) (Z:={0,±1,±2,...}) are considered. We consider two classes dissipative operators with separated boundary conditions both at -∞ and ∞. For each of these cases we establish a self...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- Numerische Mathematik
دوره 124 شماره
صفحات -
تاریخ انتشار 2013